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Abstract 



The gauge-mediated SUSY-breaking (GMSB) model needs entropy production 
at a relatively low temperature in the thermal history of the Universe for the un- 
wanted relics to be diluted. This requires a mechanism for the baryogenesis after 
the entropy production, and the Affleck and Dine (AD) mechanism is a promising 
Q \ candidate for it. The AD baryogenesis in the GMSB model predicts the existence 

■ of the baryonic Q ball, that is the B ball, and this may work as the dark matter in 

the Universe. In this article, we discuss the stability of the B ball in th presence of 
baryon- number violating interactions. We find that the evaporation rate increases 



}^ • monotonically with the B-ball charge because the large field value inside the B ball 



enhances the effect of the baryon-number-violating operators. While there are some 
difficulties to evaluate the evaporation rate of the B ball, we derive the evaporation 
time (lifetime) of the B ball for the mass-to-charge ratio ~ 100 MeV. The life- 
time of the B ball and the distortion of the cosmic ray positron flux and the cosmic 
background radiation from the B ball evaporation give constraints on the baryon 
number of the B ball and the interaction, if the B ball is the dark matter. We also 
discuss some unresolved properties of the B ball. 



1 Introduction 



The gauge-mediated SUSY-breaking (GMSB) model [|I| is one of the complete models 
to solve the FCNC problem in the supersymmetric (SUSY) extension of the Standard 
Model (SM). In a typical GMSB model 0, the dynamical SUSY-breaking scale is of the 
order of lO'^ GeV, so that the SUSY-breaking masses of the order of the weak scale (m) 
are generated in the SUSY SM. This predicts a gravitino with a mass ^3/2 = 100 keV, 
which is beyond the closure limit of the Universe, m3/2 < 2h'^ keV ||^. It is very unlikely 
that such a small gravitino mass can be generated in the GMSB model, while there is an 
exceptional extension for it [Q. Also, the string moduli may supply another problem in 
cosmology, since it is expected to have a mass of the order of m3/2. 

This implies that there must be a substantial entropy production in such a way that 
they can be diluted in the thermal history of the Universe, we therefore have to con- 
sider a possible mechanism of baryogenesis at a relatively low temperature. A promising 
candidate for it is the Affleck and Dine (AD) mechanism 00. One of the important 
predictions of the AD baryogenesis in the GMSB model is the existence of a stable Q ball 
H^. It can be a candidate for the dark matter (DM) in the Universe. 

The Q ball, a non-topological soliton, is a coherent state of a complex scalar field 
P]. The existence and the stability come from a global U(l) quantum number. In the 
SUSY SM the Q ball is composed of squarks and/or sleptons with baryon {B) or lepton 
(L) numbers 0. In the GMSB model, the flat directions (p composed of the squarks 
and/or sleptons are lifted up by at most a logarithmic potential for the large fleld value 
due to the SUSY breaking 0. This leads to the existence of the stable Q ball. The 

3 

mass of the Q ball originated from the flat direction is proportional to mQi, not to mQ, 
because of the logarithmic potential. Here, m is the mass scale of a logarithmic potential 
(m ~ 10'^"^ GeV). Therefore, the baryonic Q ball, the B ball, may be stable against 
the decay into nucleons if the baryon number is sufficiently large (Q ~ 10^^), since the 
mass-to- charge ratio ujq can become less than 1 GeV. 

The AD baryogenesis, which is the natural candidate for baryogenesis in the GMSB 
model, as mentioned above, can generate such large B balls. In the final stage of the 
AD baryogenesis, the coherent state of the AD scalar field, which consists of the fiat 
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direction 0, becomes unstable and instabilities develop. The Q ball is formed as a result 
of the fluctuation glowing. This behaviour of the AD field has been justified by numerical 



simulations [IC], and the largest Q ball charge is proportional to the initial field value 
of the AD field. Then, the B ball dark matter is an important prediction in the GMSB 
model, assuming the AD baryogenesis. 

The B ball DM search has already given a constraint on this scenario. Since the B 
ball with larger baryon number becomes stabler, the B ball absorbs nucleons and emits 
an energy of about 1 GeV per a nucleon when the B ball collides with nucleons. This 
process is known as the KKST process, and it is similar to the monopole-catalyzed proton 
decay |TI|. From the BAKSAN, Gyrlyanda, and Kamiokande experimental results, the 
constraint on the B ball charge is derived assuming that the B ball is the DM, and it 
should be larger than 10^^ for m < 1 TeV [|l^]. This KKST process is suppressed by the 
Coulomb repulsion if the B ball has positive electric charge. The B ball has an electron 
cloud around it and behaves as a heavy atom. The interaction with matter is similar to 
the case of nuclearites and the B ball charge should be larger than 10^^~^°, depending on 
the electric charge \V2 . 

These bounds are loosened if the supergravity contribution to the flat direction poten- 



tial is included [1^. For a large B ball charge, the field value inside the B ball becomes 
large and the supergravity contribution to the flat direction potential is not negligible. 
In this case, for a larger B ball, its radius Rq does not increase and becomes a constant, 
~ 10/m3/2. The geometrical cross section ti'Rq is smaller than when the supergravity 
contribution is not included. 

In this article we derive another constraint on the B ball DM scenario. The baryon 
number may not be an exact symmetry in nature, and high-energy physics may violate 
the baryon-number conservation, such as in the grand unified theories (GUTs). In fact, 
the AD baryogenesis needs the (5 — L)-violating operators, which kick to start the rota- 
tion and to generate the baryon number. While the interaction with only lepton-number 
violation can work for the AD baryogenesis, the baryon-number-violating operators are 
needed to generate for the B ball to be generated in the AD mechanism. In this article, 
assuming the existence of the baryon-number-violating higher-dimensional operators, we 
evaluate the lifetime (the evaporation rate) of the B ball. The higher-dimensional opera- 



2 



tors, including 0, enhance the evaporation rate by the large field value inside the B ball. 
Especially, for a larger B ball, the field value becomes larger and the evaporation rate 
is significantly enhanced. As a result, the B ball does not necessarily keep the baryon 
number beyond the age of Universe. 

The final state in the evaporation of the B ball by the baryon-number- violating oper- 
ators depends on the quantum numbers of the B ball, the symmetries of the operators, 
and Uq. If the operator violates {B -\- L), but not {B — L), and the B ball does not have 
lepton numbers at all, the final states arc (e+,7r^), (z/, 7r°), and so on for LUo > 100 McV. 
If o^o ^ 100 MeV, the final state consists of only lepton and antileptons. Photons may be 
included there. Without the knowledge about the surface dynamics on the B ball, one 
cannot estimate the evaporation rate in the case a;o ^ 100 MeV. Since the pion emission 
from the B ball is not allowed kinematically, quarks emitted by the baryon-number- 
violating operators are bounded to the surface or inside of the B ball for a while, and 
they decay or annihilate to the lighter particles. In this article, we restrict ourselves to 
evaluating the evaporation rate in the case ojo ^ 100 MeV; we discuss what may happen 
for ojo < 100 MeV. 

If e"*" is in the final state, the energy is of the order of and almost monochromatic. 
This may change the cosmic ray positron flux. The existing positron flux measurements 
are above about 100 McV, which works for our case of ljq ^ 100 MeV. Also, when 7r° 
is in the final state, it may distort the cosmic background radiation. We find that both 
observations give more stringent constraints on the evaporation rate of the B ball than a 
comparison with the age of the Universe, assuming that the B ball is the DM. 

We organize this article as follows. In the next section, we introduce the Q ball follow- 
ing Coleman's argument, and review the profile and the quantum numbers of the B ball 
originated from the fiat directions of the squarks and/or sleptons in the GMSB model. 
Also, we show the profile of the Q ball configuration after including the supergravity 
contribution to the scalar potential. In Section 3, the final state of the B ball evapora- 
tion by the baryon-number- violating operators is discussed there for both ujo ^ 100 MeV 
and ujq ^ 100 MeV. In Section 4, the evaporation rate of the B ball is presented for 
ujo ^ 100 MeV and it is compared with the age of the Universe. In Section 5, we evaluate 
the fiuxes of the e'^ and 7 emitted from the evaporation of the B ball, and give constraints 
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on the evaporation rate from the observations. Section 6 is devoted to conclusions and 
discussion. 



2 Properties of the B ball in the GMSB model 

In this section we review the Q ball in order to fix our convention, and summarize the 
properties and the profiles of the B ball in the GMSB model. As Coleman pointed out p| , 
some scalar potential with a Uq(1) symmetry predicts the Q ball to be a non-topological 
soliton. The Lagrangian of the scalar field with the Uq{1) charge q is 

C = \d,<p\'-V{<j,). (1) 

This system has two conserved quantities, the charge Q and the energy E, 

Q = t fd'xq (0^(9^0) - (9^0^)0) , (2) 



E = fd'x{\dt<P\' + \dM' + V{<P)). (3) 



In this system, the non-trivial lowest-energy state with Q fixed, that is the Q ball, is 
derived by minimizing 

E^ = E + uo\Q-ij d?xq {<P\dt<P) - (9t0*)0)| . (4) 

Here, uj is the Lagrange multiplier. The time dependence on of the Q ball configuration 
is determined as 

0(x) = 0(a;)e-*«^*, (5) 

and E^ is reduced to be 

= Jd'x[{d,4>y + V^{4>))+ooQ, (6) 

where K;(0) = V{(f)) — q^uo'^ip'^. Then, the procedure to derive the Q ball configuration is 
by deriving a solution of for the equation of motion 

dr<P + -dr<P - - ~ = 0, (7) 
r 2, 36 
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Table 1: The renormalisable F and D flat directions in the SUSY SM with the baryon 
numbers. Q and L stand for the doublet quarks and leptons, and u, d, and e for singlet 
quarks and leptons. 



with the boundary condition dr4>{0) = (f){oo) = 0, and minimizing for uj. 

We now discuss the properties and the proflles of the B ball in the GMSB model. The 
large B ball originates from the flat direction consisting of squarks and/or sleptons in 
the GMSB model. A typical example of the renormalisable F and D flat directions is 
given as 



^2 



1 



1 



Here, u and d are singlet quarks, the upper and lower indices for colour and generation, 
respectively. The coefficients in the right-hand sides of Eqs. (H), 1/ v^, is for the canonical 
normalization of 0. The renormalizable F and D flat directions in the SUSY SM are 
summarized in Ref . . Here, we present only the flat directions with the baryon numbers 
in Table |l|. Here, Q and L stand for the doublet quarks and leptons, and e for singlet 
leptons. We suppressed the gauge and generation indices. The charge for the B ball 
composed of the flat direction follows the charge of this direction. This is important for 
determining the flnal states in the evaporation of the B ball. 

Which type of the B ball listed in Table |l] is the DM in the Universe as the result 
of the AD mechanism? For a flxed baryon number Q, the stable B ball should be only 
the lightest among those of the above flat directions. Even if the AD mechanism creates 
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heavier B balls with the same baryon number Q as the lightest one, they change their 
lepton numbers by emitting neutrinos and/or anti-neutrinos, and transit to the lightest 
one. However, the transition rate should be suppressed, and the transition time might 
be longer than the age of the Universe. This is because the potential barrier between 
the two states is very high due to the large field values. Therefore, we will not calculate 
the transition rate to the lightest state, however, discuss the evaporation of the B ball 
generically. 

Next, we show the profiles of the B ball in the GMSB model. In that model the fiat 
direction has a logarithmic potential due to the radiative correction from the messenger 
sector This implies that a larger B ball with a larger field value is energetically 
favoured for the unit charge, which leads to a stable B ball. In this article the scalar 
potential of the fiat direction (p in the GMSB model is simplified as 

Here, m is the SUSY breaking scale in the SUSY SM (m ~ 0{W-^) GeV), and M is for 
the messenger quark mass in the GMSB model. In a typical model, the messenger quark 
mass is of the order of lO^GeV. While this double-log potential realizes the behaviour 
of the exact potential for ^ M, derived in Ref. 0, it has a wrong behaviour for 
(j) <^ M . However, the Q ball properties are determined by the behaviour of the potential 
for ^ M, and this potential thus works for our purposes.^ The profile of the large B 
ball derived by this potential is well approximated to be 

(j) = 00 , (10) 

qujr 

for r < Rq{= 7i/quj), and = 0, for r > Rq. The parameters in Eq. (p!oD and the B ball 
mass determined from this profile are given by the B ball charge Q as 



quj = (27rr/o)MmM)2(Q/g)- 
^0 = {^)HmM)HQ/q)i, 



^ In Ref. they adopt a single-log potential and take M = m in the potential so that the behaviour 
of the potential for <C M is realized. We find that the field value inside the Q ball in our potential with 
M = TO is a few time larger than that in the single- log potential for a fixed Q. 
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Rq = I^^J (mM)-^(g/g)3, 

ruQ = ^{2nr]o)kmM)^Q/q)i. (11) 

In the following the flat directions are deflned so that q, uo and ujq are positive. The 
parameter ?7o is fltted as 

TfX 

r/o -4.8 log — + 7.4. (12) 
quj 

from the numerical calculation. For r ^ Rq, where (j) ^ mM /qu, the above approximated 
solution is not valid. However, such a region does not contribute dominantly to the 
properties of the B ball. In the following we call the B ball given by Eqs. ([10[) and (p+f) 
as the GMSB B ball. 

From Eqs. (0), the energy per unit charge ci;o(= mg/Q) is 

loq = ^uj. (13) 

The physical meaning of u is the energy for increasing or decreasing the B ball charge by 
the unit charge, since (mg — mq^i) = uj. 



The fleld value 0o is proportional to Q*, as in Eq. ([Tl|) , and for larger Q the super- 
gravity contribution to the scalar potential may be important [Q. The scalar potential 
from the supergravity contribution is 

^SUGRA = m2/2|0|2(^l + inogl^^, (14) 

where 7713/2 is the gravitino mass and Mq is the reduced Planck mass. Here we assume 
the minimal supergravity for simplicity, and this potential becomes dominant when is 
larger than mM /m^i2- The second term in the bracket of the right-hand side comes from 
one-loop correction. If it comes from the gauge interaction, K is negative and 0(10^^). 
The existence of the Q ball solution requires negative K [|I5]. In a limit where Vgmsb is 
negligible, the Q ball conflguration is exactly given by a Gaussian form: 

= 0oe (15) 
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This solution leads to 



qU = 7713/2, 



00 = (27r5)~2|is:|3m3/2(Q/g)2, 
mq = m3/2(Q/g), (16) 



up to 0{K) 0, and then 



uo = uj. (17) 



We refer to the B ball given by Eqs. (|1^) and ([Tq) as the supergravity B ball. 

In Fig. |I| we show the mass-to-charge ratio of the B ball, ujq, as a function of Q. The 
solid line is for the B ball that comes from the flat-direction potential of the GMSB model 
(the GMSB B ball). The dashed lines are for those from the supergravity potential (the 
supergravity B ball). Here, we take m = 1 TeV, M = 10^ TeV, and q = 1/3, assuming 
that the B ball comes from udd. If taking larger m or M, the solid line is shifted to 
the upper side. For smaller q, both the solid and dashed lines go up. In the region 
Uq > 1 GeV, the B ball is unstable and decays into nucleons. While the result of the 
direct search for the B ball DM depends on the electric charge, the region for Q < 10^2 IS 
excluded for any electric charge [|l^ . 



Before finishing this section, we comment on the electric charge of the B ball. This 



should be below the maximal electric charge [16| and negligible with respect to the baryon 



number. However, there is no symmetry to forbid the B ball to have an electric charge. 
For example, in a flat direction of Eq. (^), zero electric charge means that the field values 
of Ml, di, and (J2 are equal to each other. However, this equality might be violated by 
SUSY breaking and the B ball may have a non-zero electric charge. In order to determine 
it, we also need to know the details of the scalar potential. 
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3 Evaporation of the B ball by the baryon-number- 
violating operators 



In next section we will evaluate the evaporation rate of the B ball by the baryon-number- 
violating operators, using a technique described by Cohen et al. |]T^. However, this 
technique is not applicable to cases including scattering, annihilation, or decays of the 
fermions bounded inside the Q ball. The energy release in the evaporation per unit charge 
is uj. This means that we cannot evaluate the evaporation rate for uj ^ 100 MeV, where 
no mesons can be emitted from the B ball. Here, we discuss what may happen for both 
cases, UJ ^ 100 MeV and uj ^ 100 MeV, while we will evaluate the evaporation rate only 
for UJ ^ 100 MeV in the next section. 

We assume the SUSY SM with the R parity conservation. Then the baryon-number- 
violating operators are given by F terms of the effective operator with dimension larger 
than 5 or D terms with dimension larger than 6. The lowest baryon- number- violating 
operators in the F terms are given as 

C = -^QiQjQkLi\g2 + ^UiCjUkdile^ +i^-c., (18) 

where i,j,k,l are for the generations {i ^ k). These interactions change the charges of 
the B ball by A{B - L) = and A{B + L) = -2. 

Let us discuss typical examples. First, we assume that the B ball is composed of the 
flat direction in Eq. (^, U2did2. Inserting Eq. (|^) in the above operators, we get the 
interaction^ 

3^0*2^161. (19) 

Then, the final state is (7r~e"'") for one baryon number. Here, di in the pion is supplied 
by the surface of the B ball. The production (vr^e^) is allowed from the kinematics if 
C(j ^ 210 MeV, since the typical momentum for each parton is ~ \^uj. (This will be shown 
in the next section.) 

Here, notice there are other interactions, such as 

3^0*'w3e3. (20) 
^ Notice that we define that the charge of is positive. 
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Because of kinematics, the primary fermions in this interaction, and 63", cannot be 
emitted from the B ball. The primary fermions may be bounded on the surface or inside 
of the B ball for a while, and decay or annihilate into the lighter states. In fact, U3 and 
63 behave as a massless particle inside the Q ball if the above operator (EDf) is negligible. 



In the conventional SUSY-GUT, the baryon-number-violating dimension-5 operators are 
proportional to the fermion masses [|l8l; this process may then be enhanced and dominate 
over the others. If this process dominates, the final state and the spectrum may be 
different from the one mentioned above. We need a technique to calculate the transition 
rate of the primary quark, bounded in the B ball, to lighter states in order to derive the 
evaporation rate of the B ball. Keeping this possibility in mind, we will continue our 
discussion. 

If < 210 MeV, all quarks are bounded inside or on the surface of the B ball. This sit- 
uation is also similar to the case mentioned above, and we cannot evaluate the evaporation 
rate. The final states are expected to be (e"*", e~, u, z/, u) or (7's, u, z/, u). 

Next, if the B ball is composed of the fiat direction uieiU2di, the final states are (vr*^, 
e"*", e^) or (2 tt*^) for lu^ 280 MeV. On the other hand, (vr^, z/, v) is also included in the 
final states for the B ball composed of Q1Q1Q2L1. The typical momentum of each fermion 
is uj/ 4 (the momentum of the pion is double of this). These arguments are applicable to 
the other B balls. The exceptions are cases where the B ball consists of 

uuuddde, uudddd. (21) 

In these cases, the dimension-5 operators are not effective, since the interaction with 
AB = ±2 is needed, and the dimension-7 operators with AB = 2 may work well, 

^ = -^UiUjdkdid2d3, (22) 

where i 7^ j. The final states are (27r°, e~), (vr^, vr^, e~) for the former in Eq. (pT]), and 
(27r°), (Tr"*", 7r~) for the latter if the kinematics allows the processes. 

Finally, we comment on the baryon-number-violating operators in the D terms. The 
lowest ones are 

^ = ■^d^'^^QMe^P + ^^^^^QQle^P- (23) 
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Since the relevant terms in the D-teim operators have a derivative of the scalar, the 
evaporation rate is suppressed by uj/Mq, with respect to the F-term operators. 



4 Evaporation rate of the B ball 

In this section we present the evaporation rate of the B ball obtained by the baryon- 
number-violating operators. The technique to calculate it was developed by Cohen et 



al. • They evaluate the evaporation rate of the L ball to neutrinos by the lepton- 
number-conserving interaction. The evaporation process to neutrinos is equivalent to the 
neutrino pair production on the L ball background. They construct a quantum field theory 
preserving a symmetry on the L ball background, a simultaneous time translation and L 
phase rotation, and derive the evaporation rate through the Bogoliubov transformation 
between the creation and annihilation operators in the asymptotic fields of neutrino at 
t ±00. The result is given by the transition amplitude to the outgoing anti-neutrino 
from the incoming neutrino. 

We generalize their result and apply it to the B ball evaporation. In our case the 
interactions to create the fermion pairs are baryon-number-violating. However, we may 
use the formula as the zeroth order of the Yukawa coupling constants. Also, if the inter- 
actions are the dimension-5 operators, they preserve {B — L). Then, when the B ball is 
composed of udd, we can use their result by regarding the B ball as the (B— L) ball. In 
this section, we first show the properties of the Q ball evaporation process. After that, we 
present the evaporation rate of the B ball and compare it with the age of the Universe. 

Since the generalization of the technique developed by Cohen is straight-forward, we 
summarize the result without repeating their calculation here.Q The Yukawa interaction 
contributing to the evaporation is 

C = HiiJ2 (24) 

where is replaced by the Q ball background. The global Uq(1) charge for the scalar 
is 1, while those for fermions, tpi and ?/'2, are p and {—p — 1), respectively. This Uq(1) 
symmetry stabilizes of the Q ball. When adopting the thin-wall approximation for the Q 



Eq. (3.4) in Ref. has a typo. The second term in the bracket has wrong signs. 
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ball background, is taken as 

f V e"*"* (r < R 



(r > Ro). ^^^^ 



In this set-up, the evaporation rate of the Q ball is given to be 

dQ ^ f-^ dk 



J=l/2 



(26) 



where k and j are the energy and the total orbit momentum for the fermion. The explicit 
form of the transition amplitude T{k,j) is given to be 

Wo - k - h^llj^-)ih^^lj-+ - h^'lu+)w'_ 



Tik,j] 



The transition amplitude is calculated assuming that ipi and ip2 are massless. Here, 

1,2), 
(* = 1,2), 



(27) 





= h^UikRQ) ( 


ht 


= h^Ui(u^-k)RQ) 


j±+ 


= Jj±^{k+RQ), 


j±- 


= Jj±^ik_RQ), 


< 


= k + k± — u, 



(28) 

with kj- = uj/2 ± \Jik — uj/2Y — v"^. The functions, h^j\x) {i = 1,2) and jj{x), are the 
Hankel and Bessel functions. The amplitude T{k,j) has following properties: 

i) T{0,j) = T{u,j) = 0, 

ii) T{k,j) = T{uj-k,j), 

iii) T{k,j) is independent of the fermion charge p. 

Also, if Rq is not so large with respect to cu"^, the contributions from the higher j modes 
dump quickly, and the numerical calculation is not so difficult .0 
^ In Ref. the evaporation rate of the Q baU with finite Rq is calculated. 
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In Fig. 1^ we show the evaporation rate of the Q ball, dQ/dt, for Rquj = 1 and 10 as a 
function of v/u. Here, the evaporation rate is normalized by the maximum value, 

fdQ\ u^Rl 



(29) 



The qualitative behaviour of the evaporation rate is easy to understand. When v/u > 1, 
the fermion pairs are produced only on the surface of the Q ball, since the production 
inside the Q ball is suppressed by the Pauli blocking. Then, the evaporation rate is 
bounded by the maximum outgoing flux of the fermion pair with the total energy u, and 
the maximum evaporation rate is derived as in Eq. (pQ]) |[T7[ . 



On the other hand, if v/u < 1, the fermion pair production from the outer shell of 
the Q ball, whose width is ~ 1/v, contributes to the evaporation. This is because the 
penetration length of the fermion inside the Q ball, ~ 1/v, is larger than l/u. Roughly 
speaking, regarding the decay rate of a quantum with unit charge u, the evaporation rate 
is ~ X (a;t>^) x (Rq/v) = vu'^Rq, and it is suppressed by v/u compared with {dQ/dt)Max- 
Here, uv"^ is the charge density inside the Q ball. This can be proved explicitly in the 
limit of a large Rq In this limit, the evaporation rate becomes 

* = (30) 

dt 16 ^ ' 

This behaviour is not valid for 1/f ^ Rq. In this case, the whole region inside the Q ball 
contributes to the evaporation, and then the evaporation rate behaves, as ~ v'^u'^Rq as 
in Fig. |. 

The above qualitative behaviour can be seen in the energy spectrum of the fermion in 
the evaporation. In Fig. |^ we show the fermion energy spectrum for v/uo = 10, 10~^, 10~^, 
10~^, the in cases Rquj = 1 and 10. The spectra for v/u = 10 are almost independent of 
Rquj. On the other hand, the spectrum has a peak around uj/2, and it becomes steeper 
for smaller v and larger Rq. 

Now we have prepared for calculating the evaporation rate of the B ball. Here, we use 
the parametrizations given in Eqs. (0) and ([T6|) for the thin-wall approximation. This 
approximation may make an 0(1) error for the evaluation of the evaporation rate. 

In order to make our discussion clear, we assume that the B ball is composed of the flat 
direction U2did2 (Eq. (S)), and that evaporates into (vr^e^) by the interaction in Eq. 
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In this case, the 0, which is canonically normahzed, has the baryon number q = 1/3, and 
Eq. (0) becomes 

C = ^e'i"*Miei (31) 
3M^ ^ ' 

inside the B ball. We can apply the above formula by inserting 



(32) 



and rescaling u due to the charge definition. As mentioned in the last section, ef and 
Ml in the final state share 2uj/3 from this interaction. The di in vr" is supplied from 
the surface, since di is heavy with the gaugino inside the B ball. The di shares 2uj/3 
with the other di in the final state. The time scale of the di emission is of the order of 
3ci;~^/2, and the effect on the evaporation rate is negligible.^ Here we have to note that 
TT^ is a pseudoscalar, and a chirality flip is required to create a vr" from di and Ui. The 
associated suppression F may be of the order of {mg/ f.„y ~ 10"^, where is the pion 
decay constant, because it comes from the pion current interaction J^9^vr//,r- In this 
article we do not attempt to estimate this suppression factor precisely. We take F = 1 
when we present our numerical result. 

In Fig. ^ we show the evaporation time [dQ/Qdt)^^ of the B ball composed of U2did2- 
The two solid lines are for the cases the B ball originates from the GMSB or supergravity 
scalar potentials. Here, for the supergravity B ball, we flx 7/13/2 = 300 MeV and l-ft'l"^ = 
10 and use Eq. (|16D for 0, Rq, and u. Also, for the GMSB B ball, we use Eq. (|Tl]) with 
m = 1 TeV and M = 10^ TeV. We plot the line of the GMSB B ball for u > 210 MeV. 
The suppression factor of the dimension-5 operator is taken to be 10^*^ GeV, which is 
completely beyond the constraint from the negative search in proton decay pO|| . 

In this figure, the evaporation time of the supergravity B ball decreases as Q increases 

for 10^^. This is because v is much smaller than u and Rq^, and v is proportional 

to Q. On the other hand, when 10^*^, the evaporation time increases since (dQ/dt) 

is independent of v for v/u > 1. If the evaporation time is less than the age of the 

Universe (to — lO^^years), the B ball cannot survive to this day. In this figure we take 

^ The emission of di is similar to the case oi v/u > 1 in Fig. ^; the emission rate is understood from 
the analogy. 
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^3/2 = 300 MeV. For smaller m3/2, beyond the validity of our calculation, if the qualita- 
tive nature does not change drastically, {dQ/Qdt)~^ may become larger by l/m|^2 when 
(5< 10^^(7713/2/1 GeV)"^ and l/ms/s when Q > 10^^(7713/2/1 GeV)"\ The evaporation 
time of the GMSB B ball in this figure is proportional to Q^^. This is because v <^ lu 
and V is proportional to Q^, not Q as the supergravity B ball. 

In Fig. H we show the evaporation time of the B ball composed of U2did2 as a function 
of Q and Mg, in a case where the B ball originates from the supergravity scalar potential, 
with m3/2 = 300 MeV and \K\-^ = 10. The thin lines are for (dQ/Qdt)-^ = 10°, 10^°, 
10'^'' years. The bold lines are for the age of the Universe (10^° years) and the big-bang 
nucleosynthesis time (1 second). For 10^^~^°, 0o is larger than the Planck mass, and 
the region may be disfavoured from the theoretical point of view. The region for g < 10^3 
is an unphysical region, since the B ball should be an unstable GMSB B ball from Fig. |l|. 
This means that should be larger than 10^^ GeV at least so that the B ball behaves 
as the DM of the Universe when 7773/2 = 300 MeV. 

If the B ball evaporates before the time of the big-bang nucleosynthesis, it has no effect 
on the cosmology, except for the dilution of the baryon number of the Universe. On the 
other hand, if the B ball evaporation occurs after the nucleosynthesis and the abundance 
is not negligible, the energetic particles in the final state may destroy the success of the 
big-bang nucleosynthesis. Also, the entropy production may change the expansion rate of 
the Universe. In this article, we do not discuss the AD baryogenesis in detail; however, the 
evaporation of the B ball by the baryon-number-violating operators may give a constraint 
in the AD baryogenesis on the GMSB model, which predicts a B ball with large Q. 

Finally, we comment on the other B balls. If the B ball is composed of Uidid2, the 
fermions Ui, di, and ^2 are massive inside the B ball, and this might be in conflict with 
the assumption in the above formula, Eq. (^). However, the fermion associated with the 
fiat direction, which is a linear combination of Ui, di and ^2, is still massless on the fiat 
direction condensation, and the above formula is thus still applicable. 
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5 Constraint from observation of the cosmic rays 



Even if some relic particle has a lifetime longer than the age of the Universe, the de- 
cay products may distort the cosmic-ray background. This leads to a constraint on the 

are con- 



decaying DM. In fact, the relic axion |21[] and the relic Kaluza-Klein graviton p2 



strained from the the cosmic diffused gamma ray. In this section we derive constraints on 
the B ball evaporation rate from the cosmic diffused gamma background and the cosmic 
ray position flux. 

First, we start from the position flux induced by the B ball evaporation. The positrons 
in the evaporation of the B ball have an almost monochromatic energy spectrum if the 
positron is the primary fermion in the evaporation process and if Rquj^ 1 and v/uj <^ 1, 
as discussed in the previous section. However, the positrons are diffused by the galactic 
magnetic field and lose their energy through the inverse Compton and the synchrotron 
processes, by the starlight and the cosmic microwave background. In this article, we 
consider the standard diffusion model for the propagation of positrons in the galaxy. 



which was summarized in Ref. |2^. In that article, the positron flux from the neutralino 
annihilation in the Halo was calculated, assuming that the neutralino is the DM in the 
Universe. 

In their formulae, the diffusion zone of the positron in our galaxy is a slab of the 
thickness 2L ~ 6kpc, and the positron density becomes zero outside that region, since 
the positrons escape freely there. The standard diffusion-loss equation for the positron 
density spectrum {drie+fdE) is 



d drie-i 
dt de 



K{e,x)^ 



dUe^ 

de 



d 

^We 



6(e,x; 



drif. 



de 



H — 

dtde 



(33) 



where e = E/{1 GeV). Here, if (e, x) is the diffusion constant, &(e, x) the positron energy 
loss rate, and {dnf+/dtde) the source term. 

In the diffusion zone the diffusion constant K{e) is 



Kie) = Ko{C + e' 



3x10^^ (3°-6 + e' 



0.6^ ^^2 g-l 



(34) 
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for £■ ^ 3 GeV, and the positron energy loss rate 6(e) is 

b{e) = ^e' 

~ 10~Vs-\ (35) 

By deriving the stable solution for the diffusion- loss equation in the above environment, 
the positron spectrum originated from the DM is given as 



Here, the energy- dependent diffusion time TD{e,e') is 

f°° ~ ( Irr' \ (r-r')^ 

X dr'r'f{r')h[j^^—^\e ^'^o^e^^ (37) 

where r and z are cylindric coordinates for the position of the solar system in our galaxy 
(r = 8.5 kpc and z = 0). The function /o(x) is /o(x)e~^' with /o(a;) the modified Bessel 
function, and Erf(x) is the error function. We neglect the energy- dependent part of the 
diffusion constant, and 



Av^C[--h. (38) 



The function /(r) is 



/(r) ^ J^^dz = g^UV^^^), (39) 



with goMirsph) {.^Iph = + z^) the DM density profile. The exponent m is 1 for the 
decay process of the DM, and 2 for the annihilation process. In this article, we use the 
isothermal distribution for the DM density profiles. 



,2 



gTiM{rsph) = 2 I 2 ' (^0) 

with a = 5 kpc. The A^-body simulation suggests the cuspy density profile at the centre of 
the galaxy |2^, and it may not be consistent with the isothermal distribution. However, 
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the decay process is not sensitive to the density profile, and we thus use the isothermal 
distribution for the DM here. 

If the positron flux comes from the evaporation of the B ball by the baryon-number- 
violating dimension-5 operators, the source term {dnf+ / dtdE) is 

Here, we simplify the spectrum of the positron as a monochromatic one. This is justified 
for Rquj > 1 and v <^ u giving a long evaporation time, as shown in Fig. |^. We take the 
number density of the B ball no as QuqUq = 0.3 GeV cm~'^. The diffusion and the energy 
reduction in the propagation makes a tail below quj in the positron spectrum. In Fig. ^ 
we show the primary positron flux spectrum {d^^^+/dE), which is given as {d^^^+/dE) = 
{c/ATT){dne+/dE), assuming qu = 300 MeV, q = 1/3, and (dQ/Qdt)^^ = 10^^ years. Also, 
we take c^o = as for the supergravity B ball. If cjq = 4^7/3 as for the GMSB B ball, the 
primary flux is reduced by 3/4. 

While the cosmic ray positron flux is measured for an energy larger than about 70 MeV 



25 1 , the theoretical estimate of the background has uncertainties. The secondary positron 



flux, which comes from the nuclear interaction of the primary cosmic rays in the interstellar 

as 



space, is estimated in Ref. [^. The result is fitted in Ref. |23 
d^['} 4.5e0-7 



cm"^ s'' GeV-\ (42) 



dE 1 + 650e2-3 + 1500e4-2 

The qualitative behaviour of the positron fraction (e"'"/(e^ + e~)) in the cosmic ray, in- 
creasing of the positron flux at lower energy, is realized by assuming that the positron is 
of a secondary origin. However, the effect of the solar modulation from the solar wind 
and the magnetosphere to both the positron and electron fluxes is larger at lower energy. 
Especially, the introduction of the charge-dependent solar modulation makes the fit to the 
observation worse |^ . Therefore, we derive a conservative constraint on the evaporation 
time of the B ball by imposing a condition that the peak of the primary flux from the B 
ball be smaller than ten times of the secondary flux, Eq. (^2]), and we obtain 

f dQ\~^ -, ^19 f Q^oriQ \ f quj \ 



Qdt 



> 2 g X 10^^ years ( A ( ) , (43) 

~ ^ ^ V0.3 GeVcm-sy VlOOMeV/ ' ^ ^ 
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for quj ~ 100 MeV. From Fig. this constraint means that M5 should be larger than 
10^^ GeV for Q ^ 10^^, assuming the B ball is the DM of the Universe and ^3/2 = 
300 MeV. 

Next, let us consider the distortion of the cosmic background radiation by 7r° from 
the B ball evaporation. Although the gamma, which comes from -k^ 27, is almost 
monotonic at production time, the energy is reduced by the red-shift. As a result, the 
energy spectrum of the gamma from the B ball evaporation is 



dn^ dQ 



3 



3n„^t„£.(,.)-., (44) 

for E < qu), assuming {dQ/Qdt)~^ ^ to- Here, we assume that the B ball evaporates 
through the baryon-number- violating dimension-5 operators and that the tt" has an energy 
2quj. From this equation (|4^) , the flux of gamma rays for E < qu is 

"^^^ 7.9 X 10^ cm-2 s-i sr-i GeV-^ 



dE 



( to \f QuJono \/ Up qu yl / E y 

\{dQ/Qdt)-^) \ PC yvo.iGevy vo.iGevy VO.lGeVy ^ 

Here, pc is the critical density of the Universe {pc = 1.1 /i^ x 10^^ GeV cm~'^). In Fig. ^ 
we show the spectrum of the gamma rays flux from the B ball evaporation [d^^/dE), 
assuming qu = 250 MeV, g = 1/4, and 10^^ years for the evaporation time {dQ/Qdt)~^. 
Here we fix /i = 0.7 and to = 10^° years. Also, we take ujq = u as for the supergravity B 
ball. 

EGRET determines the extragalactic gamma ray background spectrum between 0.1 
as 

^ — /"7 Qo _i_ n Q/1^ 1 / \ 



and ~ 50 GeV \27 



dE 



(7.32 ± 0.34) X 10^^ — cm-^ s"^ sr"^ GeV"\ (46) 

^ ' V0.451GeV/ ^ ' 



Here, imposing the condition that the peak of the gamma spectrum from the B ball be 
smaller than Eq. (|46|) leads us to a constraint on the evaporation time of the B ball: 

f # V > 5 X 10" years ( -j^] ( ^] (^^) " , (47) 
\QdtJ ~ ^ ^ \^ 1010 years/ V Pc / VlOOMeV/ ' ^ ' 

for quj ~ 100 MeV. This is one order of magnitude stronger than the constraint from the 



positron flux. 
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6 Conclusions and discussion 



In this article we discuss the stabihty of the B ball in the GMSB model. The B ball is 
predicted to exist in the GMSB model, assuming AD baryogenesis. While the stability 
of the B ball comes from the baryon-number conservation, the baryon-number-violating 
interaction is required to make the baryonic AD condensation, which is a seed for the B 
ball. The B ball could therefore be unstable. We find that a larger B ball evaporates 
faster, since the field value inside the B ball enhances the evaporation rate. 

The evaluation of the evaporation rate of the B ball suffers from various difficulties, and 
we therefore restrict our calculation to a B ball of the mass-to-charge ratio uq ^ 100 MeV. 
We derive the constraints on the B ball charge and the interactions, when the B ball is 
composed of ^2(^1^2, as an example. Neglecting the chiral suppression from the final state, 
the suppression factor of the baryon-number-violating dimension-5 operator M5 should 
be larger than 10'^'^(Q/10'^^)^ GeV for 1713/2 = 300 MeV, so that the evaporation time 
{dQ/Qdt)~^ is longer than the age of the Universe (to — 10^° years). The final states 
of the evaporation may include the almost monochromatic positron or vr^. If the B ball 
is the DM of the Universe, the evaporation may give a contribution to the extragalactic 
gamma ray background spectrum and to the cosmic ray positron fiux. From the current 
data, we give constraints on the evaporation time from the primary positron fiux and the 
gamma ray background, for c^o ~ 100 MeV, of {dQ/Qdt)-^> 10^^ or lO^^ years, which 
are stronger than the constraint from the age of the Universe. 

The instability of the B ball will have some impact on the cosmology of the GMSB 
model. In this article, we do not discuss detail of the AD baryogenesis. The existence of 
the B ball may be preferred, as far as the AD field with the baryon numbers condensates 
in the GMSB model. If the B ball has an evaporation time shorter than the age of the 
Universe, the evaporation may supply the entropy production and the injection of a high- 
energy positron after the big-bang nucleosynthesis. Note that the B ball may explain the 



reionization of the Universe, indicated by the Gunn- Peterson test |28|. It is well-known 
that the Universe is ionized at z ~ 3, from the absence of the strong Layman- alpha 
scattering light. While this may be explained by the huge star formation at the time, it 
may come from the late-time decay of an exotic relic particle, and the B ball may work 
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for it. 
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Figure 1: The mass-to-charge ratio of the B ball in the GMSB model. The horizontal 
axis is for the baryon number. The solid line is for the B ball, which comes from the 
flat-direction potential of the GMSB model. The dashed lines are for those from the 
supergravity contribution to the potential (m3/2 = 10~\ 10~^, and 10~^ GeV). Here, 
we take the SUSY breaking scale in the SUSY SM m = 1 TeV and the messenger scale 
M = 10^ TeV. In the region above 1 GeV the B ball is unstable. The region for Q < 10^^ 
is excluded by the direct search for the B ball dark matter. 
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Figure 2: The Q ball evaporation rate for Rqu = 1 and 10 as a function of v/u. Here, 
the evaporation rate is normalized by the maximum value {dQ/dt)uax- (See Eq. ( |29| ) for 
the definition.) In this calculation we adopt the thin- wall approximation as in Eq. (|25|) . 
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Figure 3: The energy spectrum of the fermion in the evaporation when Rquj = 1 and 
10. Here we take v/uj = 10, 10^^, 10^'^, 10~^. In this calculation we adopt the thin-wall 
approximation as in Eq. (^Sf). 
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Figure 4: The evaporation time {dQ/Qdt) ^ of the B ball composed of U2did2- The 
two solid lines are for the cases the B ball originates from the GMSB or supergravity 
scalar potentials. We take m = 1 TeV and M = 10^ TeV for the GMSB B ball Hne and 
7713/2 = 300 MeV for the supergravity B ball line. For the other parameters for the B ball 
configuration, see the text. We plot the lines for uo > 210 MeV. The suppression factor 
of the dimension-5 operator Mg is taken to be 10^° GeV. 



m3/2=300MeV 




Figure 5: The evaporation time (dQ/Qdt) ^ of the B ball composed of U2did2 as a 
function of Q and M5, in the case where the B ball originates from the supergravity scalar 
potential with m3/2 = 300 MeV and = 10. The thin lines are for {dQ / Qdt)~^ = 10°, 

10^°, 10^*^ years, and the bold hnes are for {dQ/Qdt)~^ — 10^° years and 1 second. 
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Figure 6: The primary and secondary positron flux spectra d^^+fdE. The primary flux 
comes from the evaporation of the B ball, assuming quj = 300 MeV, q = 1/3, u = Uq, 
as for the supergravity B ball, and an evaporation time {dQ/Qdt)~^ = 10^^ years. The 
number density is fixed by QujoUq = 0.3 GeV cm~'^. The line for the secondary positron 
flux is given in Ref. [p^] . 
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Figure 7: The spectrum of the gamma flux from the B ball evaporation (d^-y/dE), 
assuming qu = 250 MeV, q = 1/4, u = uq, and {dQ/Qdt)~^ = 10^^ years. The number 
density is fixed as QuqUq = pc with h = 0.7. The result from EGRET is also shown. 
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